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We study the effects of frozen boundaries in a Monte Carlo simulation near a first order phase 
transition. Recent theoretical analysis of the dynamics of first order phase transitions has enabled 
to state the scaling laws governing the critical regime of the transition. We check these new scaling 
laws performing a Monte Carlo simulation of the 2d, 8-state spin Potts model. In particular, our 
results support a pseudo-critical (3(L) finite-size scaling of the form /3(oo) + a\/L + a^jL? , instead 

£<1 ' of /3(oo) + 6>i/L d + 9 2 /L 2d . Moreover, we obtain a latent heat, Apbc = 0.294(11), which does not 

coincide with the latent heat analytically derived for the same model if periodic boundary conditions 

{^J , are assumed, Apbc = 0.486358 . . . 

■ 

£ : 

S: 
cn : 



o 

cn 
o 

CN 
O 



43 



X 
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I. INTRODUCTION 



The introduction of computer simulation methods has been a breakthrough in the study of phase transitions in 
lattice models. The rapid increase on the computers power has enabled to analyze with great accuracy the scaling 
laws, governed by the critical exponents, and even the corrections to these scaling laws. In such analysis, finite size 
effects must be taken into account carefully^. One of these effects, the disturbance from the boundary, has usually 
been dismissed by the adoption of periodic boundary conditions (PBC). Nevertheless, in some situations the adoption 
of periodic lattices may not be adequate, either for practical or theoretical reasons. This is the case of the free 
boundary conditions used in the analysis of free surfaces, or the so called boundary fields, used in the analysis of 
wetting phenomena. In the present paper we will focus on a particular election of the boundary conditions, the so 
called fixed boundary conditions (FBC), which have been recently applied to spin rnodelaitj and gauge modelai 

Second order phase transitions exhibit universality. For this reason, all the details of the system near the phase 
transition point become irrelevant for the critical exponents. By contrast, first order phase transitions are not 
universal and, hence, all details of a simulation must be considered carefully. This includes, in particular, the choice 
Qh, of boundary conditions. What are the appropriate set of scaling laws for a first order phase transition witb-EBC? 
■ This is the question we address-|in this paper. Starting from the theoretical analysis of Borgs and Kotecky&Q and 
the diploma Thesis of MedvedL2l on the dynamics of first order transitions, we present the finite size scaling laws 
applicable to the case of FBC and we check them performing a numerical simulation of the 2-d 8-state spin Potts 
modeta with FBC. 

The paper is divided as follows. In section II a brief summary of some recent simulations where FBC have been 
adopted serves as a motivation for a detailed analysis of the finite size scaling laws which are also presented and 
discussed. Section III is devoted to a discussion of our numerical simulation and the results that we have obtained. 
In section IV we analyze our results in the light of the scaling laws presented in Sec. II. Finally, in section V we give 
some concluding remarks. 



II. FIXED BOUNDARY CONDITIONS 



A. The motivation for FBC 



Recently, the-, so called gonihedric spin modefehave been proposed as a laboratory to study discrete versions of 
string theoriesE-3. All these simulations have been performed imposing standard periodic boundary conditions on a 
three-dimensional lattice. Nevertheless, in same cases three intersecting inner planes of spins were fixed to break the 
large energy degeneracy of the hamiltoniancrH. Due to the periodicity of the boundaries, this is equivalent to fix the 
spins belonging to the six planes of the 2-d boundary of the 3-d cube formed by the spins. Since for a certain range 
of the coupling parameter, in particular for n — 0, the transition is clearly of first order, the analysis of the finite size 
effects should have been done using the scaling laws presented in this paper. We expect that the application of the 
FBC scaling laws may overcome some anomalies recently observed in the analysis of this transitionEij. 
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TABLE I: Scaling laws for Periodic and Fixed Boundary Conditions. 





P.B.C. 


F.B.C. 


Cmax = 
Xmax(L) = 
Bmin (-^) — 


c (oo) + |i + O(^r) 
70+72^ + 0^) 

*o + + O(^) 


&(oo) + £+0(£) 
co+c 2 L [i + 0(i d - 1 ) 
eo + e2L d + 0(L d - 1 ) 



Another situation where the knowledge of the FBC scaling laws seems to be crucial is the issue of the triviality 
of lattice QED. Indeed, it has been claimed that the formation of artificial monopole structures, which close over 
the boundaries, in a simulation of thej-4-d U(l) gauge model may be responsible for turning the phase transition of 
this model from second to first orderE£l. To avoid this problem, originated probably by an incorrect choice of the 
boundaries, it was suggested to perform the Monte Carlo simulations on a lattice with a topology of an sphere. Along 
these lines, Baig and Fortcl proposed the adoption of FBC to simulate an spherical topology. Effectively, to fix all 
the variables belonging to the 3-d border to unity is the higher dimensional equivalent of converting a 2-d plane 
square lattice to the 2-d surface of a sphere by collapsing the lines of the border to a single point. Nevertheless, 
to discriminate between a first or a second order nature for a transition, an accurate analysis of data produced is 
necessary, and, in particular, this will be only possible if one knows for certain the applicable scaling laws. 



B. The scaling laws 

Although speaking properly no critical exponents can be defined for first order phase transitions, it is usual to 
define a set of characteristic exponents, together with a set of pealing laws borrowed-from those of the^cond order 
phase transitions. The pioneering work of Privman and Fishertj, Binder and LandauQ and Challa et alu, provided a 
phenomenological understanding of the scaling for first order transitions. A more rigorous theoretical justification for 
these first order scaling laws was presented by Borgs and KoteckyBO. The formulation of its applicability to finite size 
scaling expressions in terms of the lattice size was the work of Borgs, Kotecky and Miracle-SoleM and, independently, 
of JankeEj. But in all these developments the existence of periodic boundary conditions was assumed. Recently, 
though, Borgs and Koteckyu have extended their analysis to include surface effects in addition to the standard volume 
effects which govern first order transitions. Following this work, MendevE3 has deduced the scaling laws for the spin 
Potts model in the presence of surface effects, in particular adopting boundary conditions other than the periodic 
ones. 1 — 1 

Following the general analysis of MendevEl, finite size scaling laws in terms of the lattice size for the case of fixed 
boundary conditions can easily be deduced. They are summarized in Table |, together with the standard laws for 
periodic conditions. In the rest of this paper we will check these modified scaling laws with the results of our numerical 
simulation. 

It should be noticed that the suggestion that in the case of free boundary conditions every transition is shifted by a 
l/L correction term caused by surface effects is quite old. Bindercil, for instance, reports on a series of experimental 
resultaL3 supporting this conclusion. 



III. NUMERICAL SIMULATION 

To test the scaling laws of Table |. we have performed a numerical simulation of the 2-d 8-state spin Potts model 
defined by the partition function 

Z poUs = e ~" B > 
{»•«} 

where the energy is 

E =~Y, 8 ^v fa = !,•», 8), (2) 

(ij) 

with j3 — J/kT in natural units. It is well known that this model exhibits a first order phase transition!!] and for this 
reason it has been chosen as a test model in several previous studies. 
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Fixed boundary conditions have been implemented along the lines stated by Baig and Forti. In a 2-d grid with 
points labeled by (n x ,n y ), all spins corresponding to the lattice points (l,n y ) and (n x ,l) for n x ,n y = l,...,L have 
been fixed during all the simulation at its initial values a = 1. With this precaution, the structure of the program, 
that implements PBC, assures the persistence of the frozen boundary. 

We have performed the lattice updating applying a well tested head bath algorithm. During the simulation we 
recorded time series files for the energy E and the magnetization M defined as 

M = g maxW - L* 

g - 1 

where q = 8 and rii is the number of spins in a given orientation. 

Table Q summarizes the details of the simulations that have performed from L — 70 up to L = 350. The number of 
production Monte Carlo sweeps varies from n pro d = 6 000 000 for L = 70, to 7i pro d = 32 700 000 for L = 350. Since 
we took measurements only every nfn p — 8 sweeps, the number of total measurements per run is ro mea s = fiprod/^flip- 
We left at least 20 nfli p r e thermalization sweeps before taking measurements&llOEj. To estimate the autocorrelation 
time of energy measurements r e , we have applied two different methods. First, we use the fact that r e enters the error 
estimate 6jk = y2 r e /n meas e na i ve for the mean energy < E > of n meas correlated energy measurements of variance 



eLivo = E (< E > -E 3 ) 2 /{n mcas - 1). 



The "true? jerror estimate £jk is obtained splitting the energy time-series into 50 bins, which were in their turn 
jackknivedonJ to decrease the bias in the analysis. The second way of obtaining t c is by a direct computation of the 
integrated autocorrelation time 

rint _l + Vd ilk ) ^^^ U - 1 \E l -<E>)(E t+J -<E>) 

where fc max is a suitable cut-o around 6Tg nt < fc max < 10 r" 1 *. The corresponding error in r™* is derived from the 
a priori formula ^2(2 fc max + 1 )/ n meas t 1 * a . 

In Fig. [I] we present the energy time-series for the L = 300 and /3mc = 1.34146 simulation run. The expected 
characteristic behaviour for a first order phase transition can be clearly seen. The system remains in one of the two 
coexisting phases for a long period of time. The energy histogram for the full series is also presented in this figure. 
The similar height of the two peaks confirms that the simulation was performed very near the pseudo-critical inverse 
temperature. 

It is instructive to compare the energy histograms corresponding to the adoption of fixed or periodic boundary 
conditions. To this end we have performed two different Monte Carlo runs, close to the respective pseudo-critical 
inverse temperatures, which are fij^f — 1.342027 and /3|^ c = 1.3378 for a lattice size L = 100. These simulations 
has been done using 8 000 000 production sweeps, with nfn p — 8, discarding the initial 250000 (150000) sweeps in 
the case of PBC (FBC) for the thermalization of the system. Both histograms can be seen in Fig. 0. They show 
the characteristic two-peaks structure. Nevertheless, the latent heat, i.e., the separation between the maximum of 
the two peaks, is clearly smaller for Fixed Boundary Conditions. This qualitative observation suggest that a simple 
analysis of the energy histograms of a true first order phase transition simulated with Fixed Boundary Conditions 
might be misleading. Effectively, if the lattice size is not large enough the energy histogram could show (apparently) 
a single peak and, in consequence, one can get the erroneous conclusion that the model exhibit a second order phase 
transition. Nevertheless, even with FBC the evolution of the energy histograms when the size of the system increases 
shown in Fig. |^ [L = 100), and in Figjl] [L = 300), exhibit the expected behaviour of a first order transition. This 
observation may be relevant in the interpretation of the analysis of Baig and Fortcl, where a disappearance of a two 
peaks structure was observed when FBC were imposed to the systemtll 

In addition to the qualitative analysis of the histograms, we have computed the specific heat, magnetic susceptibility 
and the Binder kurtosis parameter at nearby values of (3mc by means of standard reweighting techniquesEH. They 
are defined as 
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COS) = ^{{E*)-{E)% (4) 



2 



XW) - ^((M 2 )-(M)% (5) 
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TABLE II: Monte Carlo parameters of the simulation. L 2 is the lattice size, ntherm the number of Monte Carlo sweeps during 
thermalization, and n pro d the number of production runs. Measurements were taken every nfu p = 8 Monte Carlo sweeps for all 
the simulations. 



L Pmc "-the 



70 
84 
100 
126 
150 
200 
226 
250 
278 
300 
350 



T^prod 



"therm /"flip "prod /"flip 



1.3343 100 000 6 000 000 144 128(12) 87 2604 

1.3363 100 000 6 000 000 208 240(25) 60 1803 

1.3378 150 000 8 000 000 357 394(38) 53 1441 

1.33909 250 000 8 000 000 883 847(122) 35 567 

1.3398 400 000 10 000 000 1320 1341(215) 38 474 

1.3407 900 000 12 000 000 4664 5434(1582) 24 161 

1.34102 1 200 000 16 000 000 7287 6991(1476) 21 138 

1.341205 1 600 000 18 000 000 9072 9700(2454) 22 124 

1.34138 2 200 000 18 800 000 11743 16969(5058) 23 100 

1.34146 3 000 000 22 000 000 15429 27765(12969) 24 89 

1.34162 4 000 000 32 700 000 25632 53623(29055) 20 80 



In table III we show the extrema of the magnitudes above defined, together with their pseudo-critical inverse 
temperatures. The error bars of these quantities have been estimated splitting the time-series data into 50 bins, 
which were jackknived to decrease the bias in the analysis of reweighted data. 



IV. SCALING LAWS ANALYSIS 



Once we have the results from the numerical simulation on finite lattices, we can proceed to analyze the data 
imposing the scaling laws of Table |. 



A. Analysis of the pseudo-critical inverse temperature 



In Table |V| we present the results of fitting the pseudo-critical betas of C max , Xmax and B m i n to the ansatz 
(3 C + ai/L + a 2 /L 2 suggested by the finite-size scaling laws presented in Table |. Notice that we have performed two 
set of fits, one for the full range 84 < L < 350, and a second including only results from the lattice sizes 100 < L < 350. 
Notice that the fits are extremely good even for the initial range 84 < L <Jip0, but they improve slightly if L = 84 
is discarded. Remember that reasonable fits should have a goodness-of-fitcj, Q, above 0.05. Fig. || depicts the fit 
for /3™ ax (L) in the range 84 < L < 350. The exact critical inverse temperature for the 2d 8-state Potts model is 

f3 c {exact) = ln(l + ^f(8)) = 1.342454 . . .. Our results of Table |[V] are in perfect agreement with this value. 

We have also fitted our data to the ansatz (3P eaks = /?(oo) + 6\/L 2 + 02 /L 4 corresponding to the PBC finite-size 
scaling law. Even though the goodness-of-fit, Q, obtained does not allow to discard the fits, the infinite volume /3 c (oo) 
resulting from them do not coincide with the exactly known value, showing that this ansatz is unsuitable. I.e. for 
the ^(L) in the range 84 < L < 350, the fit produces Q = 0.10 and /3 c (oc) = 1.342063(11) and for /3* ax (L) in the 
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TABLE III: Extrema for the (finite lattice) specific heat, C ma x, the susceptibility, Xmax, and the energetic Binder parameter, 
Bmin, together with their respective pseudo-critical inverse temperatures. 



L 


B c 

Mm ax 


Cmax 


B x 

rraax 


Xmax 


/3 s 




70 


1.334469(53) 


89.26(96) 


1.334212(52) 


95.3(1.2) 


1.333966(52) 


0.660221(74) 


84 


1.336360(46) 


124.6(1.4) 


1.336215(45) 


143.8(1.7) 


1.336040(46) 


0.660441(71) 


100 


1.337705(34) 


171.3(1.9) 


1.337620(33) 


210.5(2.5) 


1.337492(33) 


0.660657(69) 


126 


1.339124(33) 


268.4(4.1) 


1.339081(33) 


355.8(5.7) 


1.339000(33) 


0.660774(94) 


150 


1.339905(25) 


391.5(7.4) 


1.339881(25) 


550(11) 


1.339821(25) 


0.66058(12) 


200 


1.340747(23) 


757(16) 


1.340739(22) 


1144(25) 


1.340704(22) 


0.66006(15) 


226 


1.341046(18) 


1006(27) 


1.341042(18) 


1554(43) 


1.341014(18) 


0.65981(19) 


250 


1.341229(15) 


1285(29) 


1.341225(15) 


2033(47) 


1.341203(15) 


0.65950(17) 


278 


1.341379(12) 


1695(42) 


1.341377(12) 


2735(69) 


1.341358(12) 


0.65900(20) 


300 


1.341493(12) 


2083(51) 


1.341491(12) 


3413(86) 


1.341475(12) 


0.65856(21) 


350 


1.3416496(92) 


3176(77) 


1.3416490(92) 5338(130) 


1.3416373(92) 0.65754(23) 



range 100 < L < 350, the results are Q = 0.21 and (3 c (oo) = 1.342079(13). 



TABLE IV: Pseudo-critical inverse temperature fits. Q is the goodness-of-fit. Recall that the exact critical inverse temperature 
for the model is f3 c (exact) = ln(l + vT 8 )) = 1.342454 . . . 



range L's 


P£**(L) =Pc + ai/L + a 2 /L 2 


ffi ax (i) =Pc + a 1 /L + a 2 /L 2 


/3max(i) = Pc + a^/L + ai/L 2 




Q /3 C ai a2 


Q /3 C ai ei2 


Q Pc ai a2 


84 - 350 
100 - 350 


0.11 1.342494(38) -0.219(15) -25.5(1.1) 
0.72 1.342423(46) -0.187(19) -28.6(1.6) 


0.13 1.342478(38) -0.208(15) -27.2(1.1) 
0.73 1.342408(46) -0.177(18) -30.3(1.6) 


0.13 1.342481(38) -0.210(15) -28.3(1.1) 
0.73 1.342412(46) -0.180(18) -31.3(1.6) 



B. Analysis of C ma x, Xmax and _B min 

The results of the fits to the specific heat and susceptibility maxima, C max and Xmax, together with the kurtosis 
minimum are summarized in Table 0. As before, we show the fits for two ranges of lattice sizes. Notice that the linear 
correction coefficients, c\ and ei, are two orders of magnitude larger than the coefficients, ci and e2, of the dominant 
contribution L 2 . This makes necessary to adjust the data to the ansatz C max (L) = cq + c\ L + C2 L 2 , and allows to 
estimate the corrections to the leading term. 

In simulations with PBC, the correction to the leading term is of the order ji/L 2 . If we fit our specific heat 
data in the range L = 126 — 350 to the ansatz C max (L) = -f + ji/L 2 + 72 L 2 , the goodness-of-fit is Q — 0.0003 
with an absurdly high value for 71 . On the other hand, if we do not allow for a correction term and fit the data to 
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Cmax(X) = 70+72 L 2 ,jthe goodness-of-fit turns out to be 0. 

The work of MedvedEl shows that the coefficient of L 2 in the finite size scaling of C max is related to the latent, J 
Afbc, via c 2 = (Afbc /3 c /2) 2 . In fact, it is the same relationship that holds for periodic boundary conditionaJt 
If we use our estimation C2 = 0.0389(29) from Table ^ and /3 C = ln(l + ^/(8)), we obtain for the latent heat 

Afbc = 0.294(11). 



(7) 



Another way of estimating the latent heat is from the direct calculation, right at the transition, of the internal 
energies per site of the ordered and disordered phases, e or( j = E OI< i/V, and edis = E^/V. Of course, the latent heat is 
just A = edis — e or d- Lee and Kosterlitz proposedEj to reweight a given energy histogram until both peaks have equal 
height. The locations of the two maxima in the histogram can be taken as finite size estimates, e a (L) and ed(L\—foi 
the infinite- volume limits at /3 C of Crni and edis- The scaling of e (L) and e<i(L) for fixed boundary conditional!] as 
well as periodic pbpundary conditionsca is e (L) = e or( j + 0(1/ L) and &&(L) = edis + 0(1/ L). 

We smootheoO our energy histograms to reduce the noise and searched for e Q (L) and e^(L). Table VI shows the 
estimations that we found. Fitting them to the ansatz e a (L) = e^d + ki/L and e^(L) = edis + k^/L, we obtained 
e or d = —1.6032(48) and edis = —1.3114(92), with goodness-of-fitE3 Q = 1 and Q = 0.9 respectively. Consequently 
another estimation for the latent heat is 



A = 0.292(10). 



(8) 



The agreemen t wi th our previous estimation could not be better: it is quite comforting. 

R.J. BaxterEJO derived an analytical expression for the latent heat of the q-state-Potts model-assuming periodic 
boundary conditions. Numerical evaluations of his expression are tabulated in WuE2l and Janketj. For q = 8, the 
latent heat for the Potts model with periodic boundary conditions is Apbc = 0.486358 . . . Obviously our estimations 
of the latent heat do not coincide with this value, but it should not be so surprising in view of Fig. |^, where it can 
be seen that, for L = 100, the distance between peaks for P.B.C. is so different from the distance between peaks for 
F.B.C. Although such differences could tend towards the same value with L — » oo, our analysis indicates that in fact 
they do not. 

Notice that, unlike the latent heat, the a.Ti a|ytira .lly known infinite-volume critical inverse temperature, (3 C = ln(l + 
y/(q)), for the g-state Potts model is derivedl!30'E3'E2l using the self-duality property of the model, which is independent 
of boundary conditions when L — > oo. Let us recall that our estimations of (3 C are consistent with j3 c = 1.342454 . . . 



TABLE V: -Fits on the extrema o/C ma x, Xmax and B m i n 



range L's 


C m ax(i) = C + Cl L + C2 L 2 


Xmax(i) = e + ci L + e 2 L 2 


B min (L) = Bo + B x /L + B 2 /L 2 




Q C Cl C 2 


Q eo ei e 2 


Q So B\ B2 


100 - 350 
126 - 350 
150 - 350 


0.012 254(25) -4.24(35) 0.0342(11) 
0.15 427(65) -6.14(75) 0.0389(29) 


0.033 766(101) -11.9(1.2) 0.0691(31) 
0.38 1262(203) -16.7(2.1) 0.0798(49) 


0.011 0.65461(40) 1.52(13) -92.1(9.2) 
0.11 0.65295(72) 2.19(28) -153(24) 



V. CONCLUSIONS 



The first order phase transition finite-size scaling laws for Fixed Boundary Condition lattices of Borgs-Kotecky- 
Medved have been presented, tested and shown to be the only ones that hold for the 2d 8-state Potts model. 

It is clear from our analysis that Monte Carlo simulations for FBC are necessarily going to be much more time 
consuming than those for PBC, since for PBC the system sets into the finite-size scaling region as j3 c (L) = (3 c (oo) + 
Q\/L , while for FBC it does it at the slower pace of C (L) = f3 c (oo) + a\/L. Besides, we have found that the latent 
heat is affected by the boundaries. 

It is a pleasure to thank W. Janke and A. Salas for very stimulating discussions. Financial support from CICYT 
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TABLE VI: Finite size estimates e {L) and ed(L). They are obtained by reweighting the energy histograms until both peaks 
have equal heights. The infinite-volume ordered and disordered energies are estimated from the ansatz e a {L) — e or d + k±/L and 
e d (L) = e d is + k 2 /L. 



L 




ed 


1UU 






126 


-1.586(10) 


-1.398(20) 


150 


-1.587(18) 


-1.398(13) 


226 


-1.5965(93) 


-1.3623(91) 


250 


-1.5944(83) 


-1.362(14) 


278 


-1.5970(39) 


-1.350(17) 


300 


-1.5960(27) 


-1.3452(98) 


350 


-1.5958(31) 


-1.337(13) 


oo 


-1.6032(48) 


-1.3114(92) 



contracts AEN98-0431, AEN99-0766, is acknowledged. 
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FIG. 1: Energy time series for L = 300 and (3mc = 1.34146 
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FIG. 2: Energy histograms for L = 100 and P^c = 1.342027 and /jf/^f = 1.3378 
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FIG. 3: Finite-size scaling analysis of the pseudo-critical /3max m the range L — 84 — 350 by means of the ansatz /3max(-^) = 
Pc + oli/ 'L + a2/ L 2 . The infinite volume critical point obtained from the fit is f3 c — 1.342478(38), with a goodness-of-fit Q = 0.13. 



